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ABSTRACT 

Context. A low- or intermediate-mass star is believed to maintain a spherical shape throughout the evolution from the main sequence 
to the Asymptotic Giant Branch (AGB) phase. However, many post-AGB objects and planetary nebulae exhibit non-spherical sym- 
metry. Several candidates have been suggested as factors that can play a role in this change of morphology, but the problem is still not 
well understood. Magnetic fields are one of these possible agents. 

Aims. We aim to detect the magnetic field and infer its properties around four AGB stars using H2O maser observations. The sam- 
ple we observed consists of the following sources: the semi-regular variable RT Vir, and the Mira variables AP Lyn, IK Tau, and 
IRC+60370. 

Methods. We observed the 6i_6— 52,3 H2O maser rotational transition in full-polarization mode to determine its linear and circular 
polarization. Based on the Zeeman effect, one can infer the properties of the magnetic field from the maser polarization analysis. 
Results. We detected a total of 238 maser features in three of the four observed sources. No masers were found toward AP Lyn. The 
observed masers are all located between 2.4 and 53.0 AU from the stars. Linear and circular polarization was found in 18 and 11 
maser features, respectively. 

Conclusions. We more than doubled the number of AGB stars in which a magnetic field has been detected from H2O maser polar- 
ization. Our results confirm the presence of fields around IK Tau, RT Vir, and IRC+60370. The strength of the field along the line of 
sight is found to be between 47 and 331 mG in the FLO maser region. Extrapolating this result to the surface of the stars, assuming a 
toroidal field (oc r" 1 ), we find magnetic fields of 0.3-6.9 G on the stellar surfaces. If, instead of a toroidal field, we assume a poloidal 
field (oc r~ 2 ), then the extrapolated magnetic field strength on the stellar surfaces are in the range between 2.2 and — 115 G Finally, if 
a dipole field (oc r~ 3 ) is assumed, the field strength on the surface of the star is found to be between 15.8 and ~1945 G. The magnetic 
energy of our sources is higher than the thermal and kinetic energy in the ILO maser region of this class of objects. This leads us to 
conclude that, indeed, magnetic fields probably play an important role in shaping the outflows of evolved stars. 
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1. Introduction fields around AGB and post-AGB stars (e.g., Am iri et all 
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2011; Perez-Sanchez et~alf I20111 iLeal-Ferreira et al.l 12012 : 



Low- and intermediate-mass stars (0.8-8 M Q ) are believed Vlemmings et al. 2012). However, the sample of low and 

y , to maintain their sphericity until the asymptotic giant branch intermediate mass evolved stars around which magnetic fields 

(AG B) phase. Even tho ugh some AGB stars are sli ghtly elliptical have been measured is still small. So far, detections of magnetic 

(e.g.,|Reid_&_Menten 20071: ICasfio-Carrizol l2010l). many plane- field from H 2 maser polarization w ere reported around two 

tary nebulae (PNe) do not present any spherical symmetry. How AGB stars only; U Her and U Ori dVlemmings et al.l l2002t 

an almost-spherical AGB star gives rise to a non-spherical PN is [2005). Also, the morphology and strength of the magnetic field 

still an open question. A companion to the star (binary system as a function of radial distance throughout the circumstellar 

or a massive planet), disk interaction, the influence of magnetic envelope is still unclear. Observations of different magnetic 

fields, or a combinatio n of these agents are candidates to ex- field tracers are needed to constrain the field dependence on 

plain this phenomenon (Balick & Frank 200l lFrank et al.l2007t the radial distance from the star and, therefore, improve future 

iNordhause t al. 2007 and references therein). MHD simulations. 

Magneto-hydrodynamic (MHD) simulations show that Different maser species can provide information about differ- 

the magnetic field can be an important agen t in shaping ent regions around these objects. While SiO masers are expected 



post-AGBs and PNe (e.g., Garcia-Segura et al. 1999, 2005; to be found within the extended atmosphere of the star (between 
iGarcia-Diaz et al.l 120081: iDennis et al.l 120091) . Moreover, the photosphere and the dust formation zone), OH masers are 
recent observations support the presence of magnetic detected much further out (~65-650AU). The H2O masers emit 
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Table 1. Low- and high-resolution correlation passes 



Label N d 



BW 

(MHz) 



Av 
(km/s) 



PolMode 



Low 128 1.0 0.104 Full (LL,RR,LR,RL) 

High 512 1.0 0.026 Dual (RR,LL) 

Correlation parameters for the low- and high-resolution correlation 
passes. Description of Cols. 1 to 5: The label of the observed data - 
low- (Low) and high- (High) resolution - (Label), the number of chan- 
nels (N dra „ s ), the bandwidth (BW), the channel width (Av), and the po- 
larization mode (PolMode). 



at an intermediate distance to the star, between the SiO and OH 
maser regions. The distance of the H2O masers from the star is 
expected to lie within a few to less than a h undred AU (e.g., 
ICohen|[T987l: [Bowers et al.l[T989tlElitzud[T99l . 

The present work aims to enlarge the number of mag- 
netic field detections around low- and intermediate-mass evolved 
stars. We imaged five sources of this class using very-long- 
baseline interferometry (VLB I), in full-polarization mode, with 
the goal of detecti ng H?Q masers around them. As a result of 
Zeeman splitting (Zeeman 1897), we can measure the mag- 
netic field signature on mase r lines by investigating the pola rized 
emission of the masers (e.g., Vlemmin gs et al.ll200ll l2006). 

Our sample is composed of the pre-PN OH231.8+4.2, 
the semi-regular variable RT Vir, and the Mira variables 
AP Lyn, IK Tau, and IRC+60370. W e pres ented the results 
of OH23 1.8+4.2 in iLeal-Ferreira et al.l (120121) . The analysis of 
the four remaining sources is presented in the present paper. 
Single-dish SiO maser observati ons in full-polariz ation mode 
have been previously reported by Herpi n et al.l (120 06) for RT Vir, 
AP Lyn, and IK Tau. Their results show a magnetic field of 
< B||[G] < 5.6 in RT Vir, 0.9 < B {l [G] < 5.6 in AP Lyn, and 
1.9 < B||[G] < 6.0 in IK Tau. The AGB star RT Vir also shows 
strong circular polarization in single dish OH maser observa- 
tions, indicating a strong global magnetic field dSzvmczak et al.l 
2001). We did not find any literature reports concerning the 
magnetic field for IRC+60370 in the SiO maser region, nor for 
AP Lyn, IK Tau, and IRC+60370 in the OH maser region. 

This paper is structured as follows: in Sect. 2, we describe 
the observations, data reduction, and calibration; in Sect. 3, we 
present the results; in Sect. 4, we discuss the results and, in 
Sect. 5, we conclude the analysis. 

2. Observations and data reduction 

We used the NRACQ Very Long Baseline Array (VLB A) to ob- 
serve the H2O 6i,6-52,3 rotational maser transition at a rest fre- 
quency toward 22.235081 GHz of the stars in our sample. In 
each observing run, we used two baseband filters and performed 
separate lower (Low) and higher (High) resolution correlation 
passes. The first was performed in full-polarization mode and 
the second in dual-polarization mode. We show the characteris- 
tics of the Low and High correlation passes in Table Q] and the 
individual observation details of each source in Table [2] 

We observed different calibrators for each target. Each cal- 
ibrator was observed during the same run as its correspond- 
ing target. For the calibration of RT Vir, we used 3C84 (band- 
pass, delay, polarization leakage, and amplitude). To calibrate 



1 The National Radio Astronomy Observatory (NRAO) is a facility of 
the National Science Foundation operated under cooperative agreement 
by Associated Universities, Inc. 



IK Tau, we used J0238+16 (bandpass, delay, and amplitude) 
and 3C84 (polarization leakage). To calibrate IRC+60370, we 
used BLLAC (bandpass, delay, polarization leakage, polariza- 
tion absolute angle, and amplitude). Unfortunately, no good ab- 
solute polarization angle calibrator were available for RT Vir and 
IK Tau, making it impossible to determine the absolute direc- 
tion of the linear polarization vectors (also referred to as elec- 
tric vector position angle; EVPA). However, the relative EVPA 
angles for individual polarized components within RT Vir are 
still correct (no linear polarization was detected for IK Tau). 
To determine the absolute EVPA of IRC+60370, we created a 
map of BLLAC and compared the direction of the measured 
EVPA with that reported in the VLA/VLBA polarization calibra- 
tion databasfl Our IRC+60370 observation was carried out be- 
tween the calibration observations of February 21 and March 19, 
2009 in that database, where the polarization angle of BLLAC 
changed from 25.7° to 26.0°. We thus adopted a reference angle 
of 25.8° to obtain the absolute EVPA. 

After an initial analysis of the raw data, we did not detect 
any maser emission around AP Lyn and so did not proceed with 
further calibration of this data set. For the other three targets, 
we used the Astronomical Image Processing Software Package 
(A IPS) and followed th e data reduction procedure documented 
bv lKemball et alj 0.995) to perform all the necessary calibration 
steps. This included using the AIPS task SPCAL to determine 
polarization leakage parameters using a strong maser feature. 

After the data were properly calibrated, we used the low- 
resolution data to create the image cubes for the Stokes param- 
eters /, Q, U, and V. The Q and U cubes were used to gener- 
ate the linear polarization intensity (P - -\JQ 2 + U 2 ) cubes and 
the EVPA cubes. The noise level measured on the emission-free 
channels of the low-resolution data cubes is between ~2 mJy and 
~6 mJy. The high-resolution data were used to create the data 
cubes of the Stokes parameters / and V, from which the circular 
polarization could be inferred. The noise level measured from 
the emission-free channels of the high-resolution data cubes is 
between ~5 mJy and ~1 1 mJy. 

The detection of the maser spots was done by using the pro- 
gram maser finder, as described by Surcis et al. (2011). We de- 
fined a maser feature to be successfully detected when maser 
spots located at similar spatial positions (within the beam size) 
survive the signal-to-noise ratio cutoff we adopted (8 cr) in at 
least three consecutive channels. The position of the maser fea- 
ture was taken to be the position of the maser spo t in the channel 
with the peak emission of the feature (see e.g., Richards et al. 
I20T1 . 

3. Results 

We found 85 maser features around IK Tau, 91 toward RT Vir, 
and 62 around IRC+60370. The maser identification and prop- 
erties are shown in Table [7] In Fig. Q]we show the spatial distri- 
bution of the maser components (depicted as circles). The size 
of the circles is proportional to the maser flux densities, and they 
are colored according to velocity. The black cross indicates the 
stellar position determined in Sect. 14.31 

Positive linear polarization detection is reported when suc- 
cessfully found in at least two consecutive channels. The linear 
polarization percentage (Pi) quoted in Table [7] is the Pi mea- 
sured in the brightest channel of the feature. The Pi error is 
given by the rms of the P spectrum on the feature spatial po- 
sition, scaled by the intensity peak. The EVPA error was de- 
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Table 2. Source observation details 



Code 


Source 


Class 


V,„-(IF1) 


V,,„-(IF2) 


Beam 


RA 


Deco 




Date 








(km/s) 


(km/s) 


(mas) 


(J2000) 


(J2000) 




(mm/dd/yy) 


BV067A* 


OH23 1.8+4.2 


pre-Planetary Nebula 


+44.0 


+26.0 


1.7x0.9 


07*42'" 16.93 1 


-14°42'50" 


.2 


03/01/09 


BV067B 


APLyn 


Mira variable 


-19.5 


-32.5 


- 


06 A 34 m 34.88 s 


+60°56'33 v 


'.2 


03/15/09 


BV067C 


IKTau 


Mira variable 


+42.5 


+29.5 


1.2x0.5 


03''53'"28.84 J 


+ ll°24'22 v 


.6 


02/20/09 


BV067D 


RT Vir 


Semi-regular variable 


+25.5 


+ 12.5 


1.2x0.9 


13 A 02 m 37.98 I 


+05°11'08" 


'.4 


03/15/09 


BV067E 


IRC+60370 


Mira variable 


^14.5 


-57.5 


0.8x0.5 


22''49'"58.88 I 


+60°17'56" 


'.7 


03/05/09 


*Presented 


in Leal-Ferreira et al. (2012) 

















From left to right: The project code (Code), the name of the source (Source), the nature of the source (Class), the velocity center position of each of 
the 2 filters (V; sr ), the PSF beam size (Beam), the center coordinates of the observations (RA and Dec ), and the starting observation date (Date). 



terminer] using the expres sion ctevpa - 0.5 crp/P x 180°/7r 
dWardle & Kro nberg 1974). The linear polarization results are 
enumerated in Cols. 8 (P L ) and 9 (EVPA) of Table In Fig. [Q 
the black vectors show the EVPA of the features in which linear 
polarization is present. The length of the vectors is proportional 
to the polarization percentage. 

To measure the circular polarization, we used the / and 
V spectra to perform the Zeeman analysis described by 
Vlemmin gs et al.1 (12002!) . In this approach, the fraction of cir- 
cular polarization, P v , is given by 



* V — ( y max vmin)/*max 

- 2 xA F - F , xBu[Gauss]/AvL.[km/s], 



(1) 



where V max and y„„„ are the maximum and minimum of the 
model fitted to the V spectrum, and I max is the peak flux of 
the emission. The variable Af-f> is the Zeeman splitting co- 
efficient. Its exact value depends on the relative contribution 
of each hyperfine component of the H2O 6i ? 6-52,3 rotational 
maser transition. We adopte d the value Ay-p = 0.0 18, which 
is the typical value found by Vlemmin gs et al.l (|2002). The pro- 
jected magnetic field strength along the line of sight is given by 
B\\ and Av^ is the full-width half- maximum of the I spec trum. 
Although the non-LTE analysis in Vlemmi ngs et al.l d2002l) has 
shown that the circular polarization spectra are not necessarily 
strictly proportional to dl/dv, using Ap.p> , determined by a non- 
LTE fit, introduces a fractional error of less than ~20% when 
using Eq. 1 . We report circular polarization detection when the 
magnetic field strength given by the model fit is > 3cr. The re- 
porte d errors are based on th e single channel rms using Eq. [1] 
(see iLeal-Ferreira et al.l l2012, Sect. 3.3, for further discussion). 
We list the P v and By results in Cols. 10 and 11 of Table 
where the positive sign on B\\ indicates that the direction of the 
magnetic field along the line of sight is away from the observer, 
while the negative sign corresponds to a direction towards the 
observer. In Figs. [2] and [3] we present the / and V spectra and 
the model fit of V spectra for those features in which we detect 
circular polarization. 



3.1. IKTau 

We observed a total of 642 H2O maser spots toward IK Tau. Of 
these, 525 spots survived the multi-channel criteria and comprise 
85 maser features around this source. In Figs. [T]l, we present 
the spatial distribution of these 85 maser features. In Fig. [T]II 
and[T]III, we zoom in on the two areas indicated in Fig.[T]I. 

We did not find linear polarization in any maser feature 
around IK Tau. However, circular polarization was detected 
in three features around this source: IK. 20, IK. 69, and IK. 84 
(see Table [7]). The magnetic field strength along the line of 



sight given by the model fits are: -147+15 mG, -96+31 mG, 
and +215+56 mG, respectively. These features are identified in 
Fig. [T] labeled according to their field strengths. 



3.2. RT Vir 

We observed 830 H2O maser spots toward RT Vir in total. Of 
these, 671 spots comprise 91 maser features around this source. 
In Fig.[T]lV, we present the spatial distribution of these 91 maser 
features and in Fig. [T]V we show an enlargement of the area 
indicated in Fig.fTJrV. 

We detected linear polarization in nine features toward 
RT Vir: RT.31, RT.34, RT.67, RT.68, RT.70, RT.72, RT.73, 
RT.75, and RT.90 (see Table [7]). Unfortunately, no good polar- 
ization calibrator was available, making it impossible to deter- 
mine the absolute direction of the polarization vectors (the rela- 
tive EVPA between components is still correct). 

The distribution of EVPA among the nine features clearly 
distinguishes two groups of masers. Six features, located within 
projected right ascension offset -40 < a ff 5 mas and decli- 
nation offset -10 < d off < 30 mas (Fig.[T]V) have EVPA be- 
tween -38° and -59°. Another group of features, located within 
25 < a ff < 35 mas and -15 < S ff < 5 mas (Fig.[T]V), also has 
a small EVPA dispersion (+38° < EVPA < +64°). 

Circular polarization was found in three features around 
RT Vir: RT.70, RT.75, and RT.90. From the fit of the V spec- 
tra, we inferred magnetic field strengths along the line of sight of 
-143+12 mG and -188+6 mG in RT.70 and RT.75, respectively. 
We note, however, that the model fit of RT.70 does not superim- 
pose the whole S-shape structure of its V spectrum. This is a con- 
sequence of the strong emission that peaks at 11.7 km/s (RT.75). 
Because of this strong emission, a higher noise is present in the 
spectra around 11.7 km/s. Therefore, we truncated the RT.70 
spectrum at velocity values lower than 12.6 km/s to minimize 
the impact of this noise on the fit. However, even with this trun- 
cation, a high noise is still present in part of the V spectrum and 
so the results from the model fit of RT.70 should be taken with 
caution. 

The shape of the V spectrum of RT.90 suggests blended 
emission. There are many free parameters to be taken into ac- 
count in fitting emission of blended features. Consequently, any 
attempt to obtain a magnetic field strength from RT.90 will not 
generate a solution that is unique or robust. However, it is impor- 
tant to emphasize that the shape of its V spectrum clearly indi- 
cates the presence of a magnetic field. As an example, we created 
a possible fit for this feature. The solution we found for this fit 
gives a magnetic field of -84 mG for the slightly more blue- 
shifted emission and +63 mG for the slightly more red-shifted 
feature. The features themselves are separated by approximately 
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Fig. 1. Maps of the maser features detected toward RT Vir, IK Tau, and IRC+60370. The size of the circles is proportional to the 
maser flux densities, and their colors show the velocity scale. The black crosses indicate the stellar positions (see Sect. 14.31 ). The 
black lines indicate the EVPA (for RT Vir, they could not be calibrated in terms of absolute EVPA), and their lengths are proportional 
to the fractional linear polarization. The magnetic field strength along the line of sight are also shown for the features in which we 
detected circular polarization. The x-axis is the projected offset on the plane of the sky in the direction of right ascension. The y-axis 
is the declination offset. The offsets are with respect to the reference maser. 
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Fig. 2. Plots of Stokes / (top; black line) and V (bottom; red line) spectra of all maser features in which we detected circular 
polarization. The blue lines show the best model fit to each V spectrum. The source, the maser identification, and the magnetic field 
strength along the line of sight given by the fit are presented in top-left corner of each plot. The x-axis shows Vlsr in km/s and the 
y-axis the intensity in Jy/beam. 5 
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Fig. 3. Spectra I and V (black curves) of RT.90. The dashed lines 
show the fit of the individual blended features. In red, we show 
a fit for I and V. 



4. Discussion 

4.1. Non-detection toward AP Lyn 

Several H 2O masers toward AP Lyn have been detected be- 
fore (e.g., Imai et al. 1997; Migenes et al. 1999; Colomeret al. 
2000; Sudou et al. 2002; Shintani et al. 2008; Kim et alj|201(" 



Shin tani et alJ (2008) monitored this and other sources from 
2003 to 2006, using the Iriki 20 m telescope of the VLBI 
Exploration of Radio Astrometry (VERA). They reported a high 
flux variation, and fit a maser light curve for Ap Lyn. The 
peak flux of the H2O masers reported in the ot her works vary 
from ~6 Jy (observed with the VLBA in 1996; Mi genes et al.l 
1999) up to ~120 Jy ( observed w ith the Kashima-Nobeyama 
Interferometer in 1992: ISudouet a l. 2002). Conservatively, we 
suggest an upper limit of 1 Jy for the flux density of our non- 
dete ction (in the raw data ) . 

Richards et al. (2012) give a detailed discussion of the pos- 
sible causes of H2O maser variability. They point out that a cor- 
relation of the infrared light curve and maser variabil ity can ex- 
ist. A lso, they disagree with previous papers (e.g., Shintani et al. 
2008) stating that no systematic relationship between maser 
brighness and the optical phase was found at the times of their 
observations. 



0.2 km/s and have widths of 0.38 and 0.4 km/s. We present this 
possible fit in Fig. [3] 

In Fig. [QV, RT.70, RT.75, and RT.90 are labeled with the 
magnetic field strength along the line of sight obtained from the 
model fits shown in Figs. [2] and [3] 

3.3. IRC+60370 

We observed 658 H 2 maser spots in IRC+60370 in total. Of 
these, 634 spots comprise 62 maser features around this source. 
In Fig.[T]VI, we present the spatial distribution of these 62 maser 
features and in Fig. [TJ VII andQ]VIII we show an enlargement of 
the areas indicated in Fig. Q] VI. 

We detected linear polarization in nine features toward 
IRC+60370. These nine features are concentrated in two small 
projected regions, with areas of ~2 mas 2 each. Six of them are 
located within -2 < a ff < mas and — 1 < 6 ff 5 3 mas 
(Fig. Q] VII; Zoom A). The other three features with linear polar- 
ization detection are located within 27 < a /f < 29 mas and 
6 < 5 off < 8 mas (Fig. [TJVIII; Zoom B). The EVPA of the 
masers in Zoom A vary from -97° to -39°, while the EVPA 
of the masers in Zoom B are between -74° and -132°. 

Circular polarization was found in five features around 
IRC+60370: IRC.44, IRC.45, IRC.47, IRC.48, and IRC.58. 
From the model fit to the V spectra, we measured magnetic field 
strengths along the line of sight of +47+3 mG, +266+30 mG, 
+331+82 mG, +273+18 mG, and -130+22 mG, respectively. 
We note that, once more, the model fit of several features does 
not superimpose the full expected S-shape of the V spectra. For 
this source three factors contributed to this: (i) the limit of the ob- 
served spectral range, (ii) features with similar spatial and spec- 
tral position, and (iii) the increase in the noise near -52 km/s, 
due to the strong feature IRC.44. Therefore, the results given 
by the fit to IRC.44, IRC.45, IRC.47, and IRC.48 should also 
be taken with caution. In Fig. [I] VII, the five features with cir- 
cular polarization detection are labeled with the magnetic field 
strength along the line of sight that we obtained from the model 
fits shown in Fig. [2] 



4.2. Spatial distribution of the masers 

The spatial distribution of the features around IK Tau shows a 
clear correlation between velocity and position. While masers 
with higher velocities (red circles) are concentrated in the west 
and southwest, the features with lower velocity (blue circles) ap- 
pear, mostly, in t he east and n ortheast (Fig.[T]I). This behavior is 
also reported bv lBains et al.l (2003). They suggest that the shell 
of IK Tau has an equatorial density enhancement. The bright- 
est masers would lie in an oblate spheroid and the plane of the 
equator would have an inclination angle i' with the line of sight 
(45° < /' < 90°). The eastern end of the polar axis would then be 
approaching us, explaining the east-west velocity segregation. 
This model also explains why the IK Tau observations made 
more than 10 years apart show a persistent east- west offset be- 
tween moderately red- and blue-shifted emission, although in- 
dividual masers do not survive for more than ~1.5 years. Our 
IK Tau data were observed i n 2009, almost 15 y ears later than 
the observations reported by Bai ns et al. (120031), and 2 4 years 
after the observations reported bv lYates & Cohenl(ll994l) . 

Bains et al. (2003) also observed a similar east-west velocity 
offset in RT Vir. Our data do not show a clear correlation be- 
tween velocity and position for this source (Fig.[TJl), but a mod- 
erate enhanced concentration of red-shifted features in the east 
is present, while the bluer features are concentrated in the cen- 
ter of the pl ot. This is differen t from the east-west relation seen 
in Fig. 6 of lBains et alJ d2003l) . In their figure, the red-shifted 
masers are located on the western side, and the blue-shifted fea- 
tures are concentrated on the eastern side. 

An individual H2O maser has its life time estimated to be less 
than 1-2 years. Multi-epoch imaging of 22 GHz H2O masers of- 
ten s hows major changes i n the maser distribution over the years 
(e.g jRichards et alJJ2012l) . IK Tau is, therefore, an exception to 
this behavior. 



4.3. Stellar Position 

Some of the analysis discussed in this paper requires informa- 
tion concerning the stellar position in relation to the observed 
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Table 3. Stellar position 





IKTau 


RTVir 


IRC+60370 




a, 6 (mas) 


a, 5 (mas) 


a, S (mas) 


i=\ 


- 


a = +04.51 


a = +07.44 




- 


6 = -13.83 


6 = +03.57 


i=2 


a = +07.43 


- 


a = +14.49 




<5 = +03.09 


- 


6 =+03.61 


i=3 


a = +28.58 


- 


tK = +11.17 




6 = -05.38 


- 


6= +01.99 


All Features 


or =+10.41 


a = +10.69 


a = +10.92 




S = +02.36 


8 = -15.32 


6 = +02.68 


Mean Position 


a = +15.47 


a = +07.60 


or = +11.01 




S = +00.02 


6 = -14.58 


6 = +02.96 



Position of the star, relative to the reference maser, for different values 
of i. The positions we obtained as the centroid of all the observed maser 
features are also shown. Finally, the mean result is reported at the bot- 
tom of the table. Columns 2 to 4 show the stellar position of IK Tau, 
RT Vir, and IRC+60370. 



masers (Sects. |4~4"1 and |4"31 >. However, the absolute stellar posi- 
tion is not known for our observations. So to infer the stellar po- 
sition , we used the shell-fitting method (lYateJl993t iBains et all 
2003). This method assumes a distribution of masers on a tridi- 
mensional sphere, with the star located in its center. All masers 
in a velocity range determined by 



V star ± i(AVus R /n) 



(2) 



are identified, where V star is the velocity of the star, AVlsr the 
total maser velocity range, and n is a number taken here to be 
equal to 8. We choose that value to restrict the selection of the 
masers to lie within a small velocity range. The constant i sets 
which ring(s) along the line of sight is considered. If i is equal to 
1 , then a ring at the same line of sight velocity as the star is taken. 
If i is bigger than 1 , then one ring in front and one behind of the 
star are considered. Once the masers are selected, the central 
position of the features is assumed to be the stellar position. We 
emphasize that the more asymmetric the maser distribution, the 
larger the uncertainty of this method. 

For each object, we varied the value of i, obtaining different 
locations for the stellar position. An additional position was cal- 
culated by taking the center point of all the observed masers. We 
assumed the stellar position to be the mean location of the differ- 
ent positions we obtained by using different values of i, and by 
using the center point of all the observed masers. In Table [3]we 
show the stellar position we calculated for each value of i and 
the mean result. 

4.4. Distance of the masers to the star 

In Fig. |4] we show, for each source, a plot of the velocity of the 
features versus their projected angular offsets from the star (6 ff, 
see Sect. 14.31 for the determination of the stellar position). For 
each source, two parabolas are fitted to the velocity-offset posi- 
tions. These fits are shown by the dotted lines in the figures. In 
the fitting process, made by eye, the area between the parabolas 
which contains all masers is minimized. The parabolas obey the 
relation 



R 



"off = 



V 



shell 



X (Vthell ~ (VLSR ~ VsmrY) 



2n1/2 



(3) 



where R is the distance to the star, V s hell the expanding velocity 
of the masers, and V star the velocity of the star. 



Assuming that the masers are located in a spherical shell 
around the star, it is possible to determine the internal (Ri) and 
external (R„) radius of this shell from the internal and exter- 
nal parabola fits, and their corresponding expansion velocities 
(v, and v ). The values we adopted for V star , the distance to the 
source, their respective references, and the fit parameters (v,-, v OJ 
Ri, R ) are shown in Table [4] 

IBains et al.l J2003I) and IRichards et ail (l201ll) also investi- 
gated the kinematics of IK Tau and RT Vir and found similar 
results for v,-, v , Rj, and R . To illustra te the co mparison with 
our results, we reproduce the fits from IBains et aD (2003 ) for 
IK Tau and RT Vir in our Fig. [4] Those authors present two 
alternative solutions for the internal fit to IK Tau. We choose 
to show only the one with the larger radius here. Their fits are 
shown in Fig.|4]by the dashed lines. We note that there is a big 
disagreement between the external fits from Ba ins et al.l (2003) 
and ours. This is probably because our observations with the 
VLBA resolve out more diffuse emissions, due to its longer base- 
lines. Additionally, our result implies that the H2O maser re- 
gions around IK Tau and RT Vir reach closer to the star than 
was determined by Bains et al. (2003). Quantitatively, we found 
Ri equal to 38 and 18 mas for I K Tau and RT Vir, r espectively. 
The fits that we reproduced from lBains et al.l (120031) correspond 
to Ri equal to 60 and 45 mas for IK Tau and RT Vir, respec- 
tively. We emphasize, however, that their alternative solution 
for the internal fit of IK Tau shows an inner radius of the H2O 
maser region closer to the s t ar than ours (Rj equal to 25 mas). 
For IK Tau, iRichards et ail d2011l) found fl, between 60 and 
75 mas for different epochs, but they also detected a faint group 
of masers with R, smaller than 64 mas (at 23 mas). For RT Vir, 
IRichards et al.l ( 20 1 1 ) found /?, between 34 and 45 mas for differ- 
ent epochs. Hence, considering the stellar radius of IK Tau and 
RT V i r to be, respectivel y, 0.8 A U and 2.8 AU dMonnier et al.l 
2004; Raglandetal. 2006; Richards et ail 1201 2|) , it seems that 
although the majority of the 22 GHz H2O masers occur outside 
a distance of -5-7 stellar radii, occasional clumps can be found 
as close as ~3 stellar radii. 



4.5. Magnetic field detection 

4.5.1 . Linear polarization and field geometry 

We measured fractional linear polarization from 0.1% to 1.4% 
around RT Vir, and between 0.2% and 1 .6% around IRC+60370. 
The non-detection of linear polarization in the features around 
IK Tau imply that, if present, it is lower than ~0.5% on the 
brightest masers. These results agree with the upper limits for 
fractional linear polarization derived from the non-detections of 
IVlemmings et all d2002l) . 

According to maser theory, the magnetic field lines can be ei- 
ther parallel or perpendicular to the EVPA. It is parallel when the 
angle 6 between the field and the direction of propagation of the 
maser is less than the Van Vleck angle (~55°), and perpendicu- 
lar wh en 8 is greater than the Van Vleck angle (Gold reich et al.l 
1973). The linear polarization is affected by 6 and the degree of 
saturation but, based on our measured values, we cannot ensure 
in which regime - parallel or perpendicular - the emission origi- 
nates. 

As shown in Sects. I3~2l and l3. 31 linear polarization has been 
detected in masers toward RT Vir and IRC+60370. In each 
of these sources, the polarized features are separated into two 
groups. In RT Vir, both of these groups show a small EVPA dis- 
persion (< 26° for all masers within a given group). If, in this 
source, we are dealing with a magnetic field perpendicular to the 
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Table 4. Distance of the masers to the star: input and output parameters 



Source 


V slar (ref) 
(km/s) 


D(ref) 

(pc) 


Vi 

(km/s) 


(km/s) 


Ri 

(mas) 


R„ 

(mas) 


Ri 

(AU) 


Ro 

(AU) 


IKTau 
RTVir 

IRC+60370 


+34.0 (K87) 
+ 18.2 (N86) 
-49.3 (108) 


265 (H97) 
133 (H97) 
1000 (108) 


4.2 
4.3 
4.0 


10.0 
10.4 
17.0 


38 
18 
5.5 


110 

135 
53 


10.1 

2.4 
5.5 


29.2 
18.0 
53.0 



From Cols. 1 to 9: the source name (Source), the velocity of the source and its reference (V slar (ref)), the distance to the source and its reference (D 
(ref)), the inner (v,) and outer (v„) expansion velo cities of the H2 O env elope, and the inner ( /?,) a nd outer (R ) dista nces o f the H2O maser r egion 
to the star, both in mas and AU. References: K87: Kirranel ( fl987T) : N86: lNvman et alj (fl986): 108: Imai et alJf2008T) : H97: lHipparcosl ( ll997l) . 



EVPA, either a poloidal or a dipole field seems to be the best 
qualitative fit of the field geometry to our results. On the other 
hand, if the field is parallel to the EVPA, the polarization vec- 
tors could trace tangent points of a toroidal field. In IRC+60370, 
the EVPA of the features have a higher dispersion, but the vec- 
tors still seem to trace a dominant direction, pointing towards the 
position of the star, especially in the features located within the 
Zoom A region. If, in this source, we are dealing with a magnetic 
field perpendicular to the EVPA, either a poloidal or a dipole 
field could be argued as probable fits to our results. On the other 
hand, if the field is parallel to the EVPA, a toroidal field may 
provide a better qualitative fit. Furthermore, we detected circular 
polarization in four features located within the Zoom A region 
and, from the model fit of their V spectra, all of them show a 
magnetic field strength with a positive sign. Inside the Zoom B 
region, however, the single feature in which we detected circu- 
lar polarization shows a magnetic field strength with a negative 
sign. These results lead to the conclusion that the component of 
the magnetic field along the line of sight points in opposite di- 
rections on either side of the star. That evidence suggests, again, 
a toroidal field around IRC+60370. 



4.5.2. Magnetic field dependence 

In Fig. [5] we show a plot of the magnetic field strength along the 
line of sight for the stars in our sample, estimated from different 
maser species, against the radial distance of these masers to the 
star. We use this plot to investigate the field dependence on R: 
B cc R~ a , where a depends on the structure of the magnetic field 
in the circumstellar envelope. When a equal to 1, it refers to a 
toroidal magnetic field, a equal to 2 corresponds to a poloidal 
field, and a equal to 3 indicates a dipole geometry. In the plot 
we show one single box where the results of OH masers occur. 
However, we emphasize that the 1665/7 MHz OH maser emis- 
sion originates in inner regions when compared to the 1612 MHz 
OH maser transition. Therefore, it is expected that magnetic field 
strength me asurements based o n the first line to be stronger than 
the second dWolak et al.ll2012l) . 

In this plot, we included polarization results of the SiO maser 
region from the literature. We took the m agnetic field strength in 
the SiO maser region from lHerpin et al.l (120061 for RT Vir (upper 
limit) and IK Tau. For IK Tau, the distance of the SiO region to 
the star was adapted from lBoboltz & Diamond! (12005). adopting 
a distance to the source of 265 pc. For RT Vir, we used a typical 
value for the radial distance of the SiO maser region (between 
2 and 5 AU from the star). Unfortunately, we did not find any 
reports of the magnetic field strength in the OH maser region 
that would allow us to make a more complete plot. For all cases, 
the major uncertainty in the plot concerns R. 

The data from RT Vir and IRC+60370 do not allow a defini- 
tive conclusion regarding the functional form of radial depen- 
dence. For IK Tau, however, even though a B oc R l dependence 



is not totally ruled out, B cc R 2 and B 
tively better fits. 



oc R 3 provide qualita- 



4.5.3. Magnetic Field on the Star Surface 

By assuming a magnetic field dependence (B oc R~ , B oc R~ 2 , 
or B oc R 3 ; see Sect. 14.5.2b . we can extrapolate the projected 
field strength to the surface of the star {B slar ). If B oc R~ a , then 



B., 



B Hl0 x (R/R*) a , 



(4) 



where Bh-,o is the field strength along the line of sight in the wa- 
ter maser region, and R t is the stellar radius (Reid et al.lll979t 
Reid 1990). However, we emphasize that the magnetic field 
could deviate from any power law if the various masers come 
from conditions with different densities, fractional ionization, 
etc. These differences in the physical conditions of the medium 
could lead to differences in how much the field is frozen in, 
dissipated, enhanced by shocks, etc. Therefore, a homogeneous 
medium is assumed in this extrapolation. 

In the analysis with Eq. [4] we investigate each source in- 
dividually, by varying the power law, with a fr om 1 to 3. For 
IK Ta u and RT Vir, we adopted R* as updated by Richards et al. 
(2012). For IRC+60370, however, there is no accurate mea- 
surement for the stellar radius. Therefore, for this source we 
adopted /?„=1.8 AU as an upper limit based on 18 //m imaging 
(Mei xneretal.ll 19991) . 

In order to define the value of R to be given as input in Eq. [4] 
for each source we created an alternative plot of velocity versus 
position offset (analogous to the procedure described in Sect. 14.41 
and Fig. 0). In these alternative plots we considered only the 
features in which we detected circular polarization, getting al- 
ternative values for /?, and R (R^ and R' ). We adopted R^ and 
R' as minimum and maximum values of R to be given in Eq. [4] 
We emphasize that these alternative plots were created with very 
few data points, and thus provide only approximate results for 
R'i and R' a . We combined R — R' t with the lowest value of Bh-,o 
that we observed (taking the error bar into account - BH 2 o mir ) to 
derive the lower limit of field strength on the surface of the star 
(B s iar„ w )- For the upper limit (B starma ), we combined R-R' with 
the highest value of B^o that we observed (taking the error bar 
into account - Bh^o,„ x )- 

In Table [5] we show the values given as input in Eq. |4](/?*, 
Bh^o,,,,,,, BH 2 o maX 7 R'-, an d R' )< an d the results ofB starmtn and B starmax 
for each source. 



4.5.4. Magnetic field energy 

One question that needs to be answered to improve our un- 
derstanding on low and intermediate mass stellar evolution is: 



Leal-Ferreira et al.: Magnetic fields around evolved stars: H2O maser polarization 



Table 5. Projected magnetic field strength on the star surface 



Source 


R, 


Bn 2 o min 


"HiOmax 


R' 


K 


"star win 


** star ma x 


"star,,,,,, 


"slar max 


^* star,„j l} 


"Starmax 




(AU) 


(mG) 


(mG) 


(AU) 


(AU) 


(G) 


(G) 


(G) 


(G) 


(G) 


(G) 














(a=l) 


(a=l) 


(o=2) 


(a=2) 


(o=3) 


(<*=3) 


IKTau 


2.8 


65 


271 


22.5 


27.8 


0.5 


2.7 


4.2 


26.8 


33.8 


-265 


RTVir 


0.8 


131 


194 


6.0 


14.6 


1.0 


3.5 


7.3 


64.9 


54.9 


-1185 


IRC+60370 


1.8 


44 


413 


12.5 


29.5 


0.3 


6.9 


2.2 


-115 


15.8 


-1945 



From Cols. 1 to 12: the source name (Source), the stellar radius (R,), the lower (B/y,o„„„) and upper (Sh,o,„ ov ) magnetic field strengths along the line 
of sight observed in the H 2 region, the input values of R on Eq.[4](/^ a.ndR' ), and the lower (B aarmjn ) and upper (6„ ,„, 0V ) limits of the projected 
magnetic field strength on the stellar surface assuming a = 1, a = 2, and a = 3. 



Table 6. Magnetic energy density 



Source 


Un 2 

(J/m 3 ) 


U Star (a=V) 
(J/m 3 ) 


U S ,ar (<*=2) 

(J/m 3 ) 


Ustar («=3) 

(J/m 3 ) 


IKTau 
RTVir 

IRC+60370 


-4.8 --3.5 
-4.2 --3.8 
-5.1 --3.2 


-3.0- -1.5 
-2.4- -1.3 
-3.4- -0.7 


-1.2-0.5 
-0.7-1.2 
-1.7-1.7 


0.7-2.4 
1.1-3.7 
0.0 - 4.2 


U (J/m 3 ) 


H 2 


Star 






nKT 


—6.2 
—5.1 


-0.5 
-0.5 







In the upper part of the table, from Cols. 1 to 4: the source name 
(Source), the log of the magnetic energy density (U = B 2 /2fio) in 
the H 2 maser region (U H2 o), in the stellar surface assuming a equal 
to one (U st ar (ff=l)), in the stellar surface assuming a equal to two 
{U s tar ( Q '=2)), and in the stellar surface assuming a equal to three (U st ar 
(a=3)). In the lower part of the table, we show the log of typical values 
for kinetic and thermal energy densities (Col. 1; Energy density) in the 
H 2 maser region (Col. 2; H 2 0) and on the stellar surface (Col. 3; Star). 



if present, does the magnetic field around evolved stars have 
enough energy to drive the shaping of these objects? 

Several magnetic field detections around AG B and post- 
AGB stars have been reporte d in re c ent years (e.g., Amiri et al] 



20 1 It iPerez-Sanchez et all l201lt iLeal-Ferreira et all 1201 2 
Vlemmings et al. 2012). To check if the magnetic energy den- 
sity (U — B 2 /2yUo) is important, we compare it with typical val- 
ues of the kinetic and thermal energy density around evolved 
stars (Table |6). For the calculation of these values we assume 
V exp ~5 km/s, n w ,~10 14 cnr 3 , and 7-2500 K at the stellar pho- 
tosphere, and V exp ~8 km/s, «h,~10 8 cm 4 , and T~500 K in the 
H2O maser region. In Table [6] we show the limits of the mag- 
netic energy density in the H2O maser region that we observed 
and the magnetic energy density extrapolated to the surface of 
the star. The limits are based on the field strengths along the line 
of sight reported in Table [5] Our results show that the magnetic 
energy density is dominant in the H2O maser region. Therefore, 
the magnetic fields probably play an important role in shaping 
the different morphologies of evolved stars that are progenitors 
ofPNe. 

The dominant energy on the surface of the star is still incon- 
clusive since different conclusions can be drawn if a dependence 
of either B oc/J -1 ,/? oc/? 2 , orB oc R^ is assumed. 



5. Conclusions 

We observed four AGB stars in order to detect H2O maser 
in full polarization at VLBI resolution. We did not detect any 



maser emission toward AP Lyn. Toward IK Tau, RT Vir, and 
IRC+60370 we detected 85, 91, and 62 features, respectively. 

A structured spatial distribution of maser velocities was ob- 
ser ved toward IK Tau . This behavior has already been reported 
by iBains et all d2003l) and an equatorial density enhancement 
model was proposed. A similar signature, but less pronounced, 
was observed toward RT Vi r, but with opposite velocity-position 
pattern to those reported bv lBains et al.l(l2003l) . 

We used the shell-fitting method to infer the projected posi- 
tion of the star relative to the observed masers. With the stellar 
position determined, we produced a plot of the angular offset of 
the masers relative to the stellar position versus the maser ve- 
locities. We fitted parabolas in these plots to determine the ac- 
tual distance of the H2O maser regions from the central stars. We 
concluded that the H2O masers we observed are located between 
10.1 and 29.2 AU from IK Tau, 2.4 and 18.0 AU from RT Vir, 
and 5.5 and 53.0 AU from IRC+60370. 

Linear polarization was observed in 18 features, nine around 
RT Vir and nine around IRC+60370. Circular polarization was 
found in 1 1 features, three around IK Tau, three around RT Vir, 
and five around IRC+60370. From a model fit of the Stokes V 
spectra of the features with statistically significant circular polar- 
ization detection, we estimated the magnetic field strength along 
the line of sight needed to generate the observed S-shape pro- 
file. The resulting projected magnetic field strengths lie between 
47+3 mG and 331+82 mG. With our polarization results, we 
more than doubled the number of AGB stars around which the 
magnetic field has been detected in the H2O maser region. 

Combining our results with published results for the mag- 
netic field measurements in the SiO maser regions, it is not 
yet possible to determine the magnetic field dependence on the 
radial distance R to the star. For IK Tau, either a dependence 
B oc R 2 or B oc R 3 seems qualitatively more likely, but B ccR 1 
is not ruled out. The results we found in the literature for RT Vir 
and IRC+60370 are not sufficient to draw stronger conclusions. 

The results we obtained for the magnetic field strength along 
the line of sight were extrapolated to the stellar surface of the 
observed sources, assuming B oc R~ l , B oc R 2 , and B oc R 3 . In 
the first case, the projected field strength on the AGB star surface 
(B star ) should be between 0.5 G and 2.7 G for IK Tau, 1.0 G and 
3.5 G for RT Vir, and 0.3 G and 6.9 G for IRC+60370. If B oc R 2 
is assumed, then B star was extrapolated to be between 4.2 G and 
26.8 G for IK Tau, 7.3 G and 64.9 G for RT Vir, and 2.2 G and 
-1 15 G for IRC+60370. If B oc R 3 , then B star was found to be 
between 33.8 G and -265 G for IK Tau, 54.9 G and - 1 185 G for 
RT Vir, and 15.8 G and -1945 G for IRC+60370. 

Finally, we compared the magnetic energy density we ob- 
served with the typical thermal and kinematic energy density 
around evolved stars. While the dominant energy density on the 
surface of the star is still inconclusive, we conclude that, in the 
H2O maser region, the magnetic energy density dominates the 
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Fig. 5. Magnetic field strength along the line of sight versus 
the radial distance of the masers to the star. The black boxes 
show typical regions of the plot where results from the litera- 
ture for SiO, H2O, and OH maser oc cur, and th ey are normal- 
ized for R,=l ( Vlemmi ngs et al.l2 002. 2005; He rpin et al.ll2006t 
Rudnitski et al. 2010). Our measurements are shown by the hol- 
low blue circles (IK Tau), hollow green squares (RT Vir), and 
filled red triangles (IRC+60370). The short-dashed, solid, and 
long-dashed inclined lines show a dependence oc R l , oc R 2 7 and 
oc R 3 for the magnetic field, respectively. The position of the 
AGB surface of a star with radius of 1 AU is also shown. 



thermal and kinematic energy density. Therefore, the magnetic 
fields cannot be ignored as one of the important agents in shap- 
ing planetary nebulae. 
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Fig. 4. Velocity of the features versus their projected offsets from 
the star. The dotted lines represent our fits; the dashed lines were 
reproduced from Bains et al. (2003). From top to bottom, we 
show the plots for IK Tau, RT Vir, and IRC+60370. 
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Table 7. Properties of the observed masers. Respectively, from Cols. 1 to 10: the source in which the maser was observed (Source), the maser 
identification (feature), projected offset on the plane of sky in direction of right ascension («<,//), offset in declination (<5<,//)> peak intensity (Peak 
Int), intensity flux (Int Flux), velocity of the peak (V pea k), Py, magnetic field strength (B cosQ), and if linear polarization was detected or not 
(LinPol). a ff and S ff are with respect to the reference feature. 



Source 



Feature 



a„ ff 

(mas) 



(mas) 



Peak Int 
(Jy/Beam) 



Int Flux 

(Jy) 



» peak 

(km/s) 



Pl 
(%) 



EVPA 
C) 



Py 
(xlO" 3 



By 
(mG) 



APLyn 



IKTau 



IK.01 


7.5 


6.6 


0.20 


0.35 


42.6 


IK.02 


3.2 


9.3 


0.08 


0.15 


42.5 


IK.03 


52.3 


-37.0 


0.06 


0.18 


42.2 


IK.04 


52.6 


-38.6 


0.33 


0.54 


42.1 


IK.05 


3.6 


-6.6 


0.20 


0.41 


42.1 


IK.06 


4.7 


-6.4 


0.17 


0.54 


42.1 


IK.07 


7.1 


3.9 


0.04 


0.09 


42.1 


IK.08 


1.3 


6.6 


0.03 


0.05 


42.1 


IK.09 


3.0 


-6.6 


0.24 


0.52 


42.0 


IK.10 


53.7 


-38.3 


0.30 


0.53 


41.9 


IK. 11 


1.8 


-4.5 


0.08 


0.13 


41.9 


IK.12 


-3.1 


-4.4 


0.22 


0.50 


41.9 


IK.13 


55.0 


-38.5 


0.27 


0.87 


41.8 


IK.14 


51.3 


^10.6 


0.21 


0.37 


41.7 


IK.15 


31.7 


-41.4 


0.14 


0.25 


41.6 


IK.16 


53.8 


^13.2 


0.12 


0.20 


41.4 


IK.17 


50.5 


-48.1 


0.18 


0.30 


41.0 


IK.18 


53.3 


-44.5 


0.07 


0.10 


40.7 


IK.19 


61.1 


-26.4 


0.15 


0.31 


40.7 


IK.20 


-0.0 


0.0 


4.93 


10.22 


40.5 


IK.21 


63.3 


-31.1 


0.04 


0.06 


40.2 


IK.22 


66.0 


-32.1 


0.07 


0.17 


40.0 


IK.23 


66.4 


-32.9 


0.06 


0.14 


39.9 


IK.24 


67.2 


-34.1 


0.04 


0.07 


39.6 


IK.25 


-16.1 


22.8 


0.06 


0.14 


38.9 


IK.26 


0.5 


9.5 


0.21 


0.42 


38.7 


IK.27 


-0.4 


9.8 


0.11 


0.26 


38.7 


IK.28 


73.5 


-16.0 


0.11 


0.28 


38.6 


IK.29 


73.0 


-17.0 


0.17 


0.92 


38.4 


IK.30 


-35.8 


45.1 


0.61 


1.06 


31.7 


IK.31 


-34.5 


41.9 


0.34 


0.77 


31.6 


IK.32 


-36.8 


46.9 


0.75 


1.44 


31.4 


IK.33 


13.9 


^10.6 


0.07 


0.18 


31.4 


IK.34 


-33.3 


40.3 


1.33 


2.54 


31.0 


IK.35 


13.9 


^12.5 


0.06 


0.13 


31.0 


IK.36 


-5.0 


13.5 


0.06 


0.10 


30.7 


IK.37 


-6.9 


13.3 


0.19 


0.59 


30.3 


IK.38 


-34.1 


44.7 


0.08 


0.15 


30.3 


IK.39 


-7.9 


12.9 


0.12 


0.45 


30.2 


IK.40 


-6.7 


7.6 


0.07 


0.15 


30.2 


IK.41 


-8.9 


12.5 


0.11 


0.24 


30.0 


IK.42 


-8.1 


8.1 


0.17 


0.57 


29.9 


IK.43 


1.0 


2.3 


0.12 


0.31 


29.9 


IK.44 


-0.8 


2.0 


1.11 


3.12 


29.8 


IK.45 


-9.2 


8.0 


0.17 


0.49 


29.8 


IK.46 


-5.8 


5.6 


0.07 


0.12 


29.5 


IK.47 


-9.0 


9.5 


0.13 


0.29 


29.5 


IK.48 


-9.8 


9.4 


0.10 


0.22 


29.5 


IK.49 


-18.0 


7.0 


0.06 


0.10 


29.4 


IK.50 


-10.5 


10.4 


0.08 


0.19 


29.4 


IK.51 


-2.9 


1.4 


0.12 


0.28 


29.3 


IK.52 


-6.7 


26.9 


0.06 


0.11 


29.3 


IK.53 


-3.7 


1.3 


0.13 


0.26 


29.2 


IK.54 


-10.6 


10.3 


0.07 


0.18 


29.2 


IK.55 


0.8 


-1.6 


0.05 


0.10 


28.8 


IK.56 


-0.0 


-2.2 


0.09 


0.51 


28.7 


IK.57 


-24.3 


5.0 


0.09 


0.20 


28.6 


IK.58 


-0.5 


-1.3 


0.15 


0.28 


28.2 


IK.59 


19.7 


59.5 


0.08 


0.15 


28.1 


IK.60 


-13.1 


3.6 


0.06 


0.11 


28.1 


IK.61 


19.1 


63.5 


0.60 


1.11 


28.0 


IK.62 


-12.0 


2.2 


0.10 


0.21 


28.0 



10.08+1.03 



-147±15 
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Source 



Feature 



a ff 
(mas) 



S off 

(mas) 



Peak Int 

(Jy/Beam) 



Int Flux 

(Jy) 



* peak 

(km/s) 



Pl 

(%) 



EVPA 

(°) 



(xlO~ 3 ) 



B||[Gohm] 

(mG) 



IKTau 



IK.63 


18.1 


62.6 


0.14 


0.21 


27.9 


- 


IK.64 


-27.6 


46.0 


0.08 


0.14 


27.8 


- 


IK.65 


-6.7 


-2.6 


0.45 


1.14 


27.7 


- 


IK.66 


-10.0 


-2.6 


0.29 


0.88 


27.7 


- 


IK.67 


-8.6 


-2.1 


0.20 


0.56 


27.7 


- 


IK.68 


-7.7 


-2.4 


0.18 


0.60 


27.6 


- 


IK.69 


16.4 


74.5 


2.97 


4.28 


27.4 


- 


IK.70 


^13.0 


-25.1 


0.81 


1.66 


27.2 


- 


IK.71 


-44.0 


-24.8 


0.32 


0.88 


27.2 


- 


IK.72 


-6.0 


-2.2 


0.16 


0.35 


27.2 


- 


IK.73 


-49.9 


-26.3 


0.23 


0.53 


27.0 


- 


IK.74 


10.0 


64.1 


0.15 


0.19 


26.8 


- 


IK.75 


19.7 


73.3 


0.17 


0.27 


26.6 


- 


IK.76 


-34.4 


-28.6 


0.12 


0.17 


26.4 


- 


IK.77 


20.9 


34.5 


0.08 


0.13 


25.9 


- 


IK.78 


55.8 


14.6 


2.37 


4.67 


25.7 


- 


IK.79 


6.9 


-51.7 


0.07 


0.15 


25.6 


- 


IK.80 


57.2 


15.1 


0.08 


0.14 


25.3 


- 


IK.81 


55.1 


15.3 


0.78 


1.96 


25.3 


- 


IK.82 


53.8 


15.6 


0.30 


0.88 


25.1 


- 


IK.83 


50.8 


15.0 


0.86 


2.50 


25.0 


- 


IK.84 


51.2 


15.7 


0.90 


2.18 


25.0 


- 


IK.85 


55.1 


15.5 


0.59 


1.29 


24.8 


- 


RT.01 


-19.3 


-1.4 


0.14 


0.16 


21.8 


- 


RT.02 


45.9 


-67.7 


0.07 


0.16 


21.8 


- 


RT.03 


47.6 


-64.9 


0.09 


0.13 


21.6 


- 


RT.04 


13.4 


8.5 


0.03 


0.04 


21.3 


- 


RT.05 


41.8 


-73.8 


0.05 


0.10 


21.2 


- 


RT.06 


-15.6 


-50.1 


0.06 


0.07 


21.1 


- 


RT.07 


8.8 


11.9 


0.06 


0.08 


21.0 


- 


RT.08 


-12.7 


-54.0 


0.03 


0.04 


21.0 


- 


RT.09 


-8.1 


-22.9 


0.29 


0.42 


20.9 


- 


RT.10 


-7.2 


-21.5 


0.22 


0.30 


20.8 


- 


RT.ll 


-28.0 


2.2 


0.04 


0.06 


20.8 


- 


RT.12 


-9.2 


-21.6 


0.12 


0.21 


20.5 


- 


RT.13 


-10.4 


-22.4 


0.26 


0.37 


20.4 


- 


RT.14 


-12.4 


-77.2 


0.04 


0.09 


20.4 


- 


RT.15 


-14.5 


-34.0 


0.64 


1.44 


20.1 


- 


RT.16 


-36.4 


-60.4 


0.07 


0.09 


20.1 


- 


RT.17 


-10.5 


-22.4 


0.21 


0.29 


19.9 


- 


RT.18 


52.6 


-72.5 


0.29 


0.40 


18.1 


- 


RT.19 


52.8 


-55.5 


0.33 


0.61 


18.1 


- 


RT.20 


49.1 


-51.8 


0.74 


1.13 


18.1 


- 


RT.21 


53.1 


-49.9 


4.65 


5.42 


18.1 


- 


RT.22 


-26.1 


-15.8 


0.25 


0.40 


18.1 


- 


RT.23 


-25.3 


6.0 


3.13 


5.03 


18.1 


- 


RT.24 


-23.5 


6.0 


0.77 


0.99 


18.0 


- 


RT.25 


53.5 


60.1 


0.18 


0.34 


18.0 


- 


RT.26 


-16.7 


-3.1 


0.90 


0.97 


17.7 


- 


RT.27 


-29.8 


27.2 


2.34 


3.73 


17.6 


- 


RT.28 


-32.4 


27.1 


0.62 


0.66 


17.5 


- 


RT.29 


-37.1 


-50.0 


0.56 


0.72 


17.1 


- 


RT.30 


-31.4 


25.6 


7.83 


9.24 


17.0 


- 


RT.31 


-35.8 


25.9 


20.96 


25.07 


17.0 


0.22±0.02 


RT.32 


-35.0 


28.0 


2.83 


3.55 


17.0 


- 


RT.33 


-9.6 


-28.8 


0.60 


0.90 


17.0 


- 


RT.34 


-17.0 


-5.4 


1.21 


1.38 


16.7 


1.41+0.04 


RT.35 


-30.8 


28.4 


0.79 


0.86 


16.6 


- 


RT.36 


-32.3 


25.1 


1.32 


3.13 


16.4 


- 


RT.37 


39.8 


-41.5 


0.32 


0.51 


16.3 


- 


RT.38 


-21.8 


10.7 


0.28 


0.35 


16.3 


- 


RT.39 


29.7 


-75.6 


0.48 


0.75 


15.9 


- 


RT.40 


90.7 


18.0 


0.68 


1.08 


15.6 


- 


RT.41 


79.5 


36.2 


0.48 


0.89 


15.6 


- 


RT.42 


86.2 


89.1 


0.44 


0.54 


15.6 


- 


RT.43 


-14.3 


-7.5 


0.59 


1.58 


15.4 


- 



5.48±1.78 



-96±31 



12.15+3.14 



+215+56 



RTVir 



-46+5 



-59+3 
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Table 7. continued. 



Source 


Feature 


"off 


Soff 


Peak Int 


Int Flux 


* peak 


Pl 


EVPA 


Pv 


B||[Gan.ij] 






(mas) 


(mas) 


(Jy/Beam) 


(Jy) 


(km/s) 


(%) 


O 


(xl0~ 3 ) 


(mG) 




RT.44 


87.8 


-2.1 


0.64 


1.21 


15.4 


- 


- 


- 


- 


RTVir 


RT.45 


-15.4 


-7.4 


1.15 


1.62 


15.3 


- 


- 


- 


- 




RT.46 


23.6 


-51.0 


0.22 


0.48 


15.2 


- 


- 


- 


- 




RT.47 


44.3 


-44.3 


0.28 


0.36 


15.2 


- 


- 


- 


- 




RT.48 


31.7 


-82.4 


0.39 


0.59 


14.9 


- 


- 


- 


- 




RT.49 


-6.7 


-8.0 


0.26 


0.56 


14.8 


- 


- 


- 


- 




RT.50 


-31.3 


27.6 


0.12 


0.12 


14.8 


- 


- 


- 


- 




RT.51 


-39.3 


20.4 


0.10 


0.13 


14.8 


- 


- 


- 


- 




RT.52 


33.9 


-90.5 


0.25 


0.36 


14.6 


- 


- 


- 


- 




RT.53 


31.8 


-36.2 


0.32 


0.35 


14.5 


- 


- 


- 


- 




RT.54 


58.5 


-16.8 


1.20 


1.33 


14.3 


- 


- 


- 


- 




RT.55 


58.2 


-50.3 


0.08 


0.11 


14.3 


- 


- 


- 


- 




RT.56 


-22.9 


-22.1 


0.08 


0.18 


14.3 


- 


- 


- 


- 




RT.57 


-20.9 


-20.9 


0.23 


0.41 


14.2 


- 


- 


- 


- 




RT.58 


2.8 


-18.3 


0.13 


0.15 


14.0 


- 


- 


- 


- 




RT.59 


29.4 


-15.8 


0.34 


0.46 


13.9 


- 


- 


- 


- 




RT.60 


3.6 


4.2 


0.58 


0.97 


13.8 


- 


- 


- 


- 




RT.61 


30.6 


-16.4 


0.15 


0.20 


13.8 


- 


- 


- 


- 




RT.62 


29.9 


-14.0 


0.35 


0.39 


13.7 


- 


- 


- 


- 




RT.63 


2.2 


1.8 


0.42 


0.48 


13.7 


- 


- 


- 


- 




RT.64 


2.4 


19.0 


0.36 


0.40 


13.7 


- 


- 


- 


- 




RT.65 


-1.9 


-1.0 


0.57 


0.71 


13.6 


- 


- 


- 


- 




RT.66 


5.2 


-54.0 


0.36 


0.49 


13.4 


- 


- 


- 


- 




RT.67 


-1.0 


-1.8 


2.08 


2.71 


13.1 


0.56±0.26 


-44+16 


- 


- 




RT.68 


0.0 


0.0 


54.53 


63.13 


12.9 


0.48+0.19 


-38+12 


- 


- 




RT.69 


8.0 


-85.1 


1.30 


1.60 


12.9 


- 


- 


- 


- 




RT.70 


30.9 


-10.6 


7.92 


8.45 


12.9 


1.12+0.28 


+49+26 


10.50+0.86 


-143±12 e 




RT.71 


-0.3 


55.2 


0.93 


1.29 


12.9 


- 


- 


- 


- 




RT.72 


-7.5 


12.0 


4.23 


5.28 


12.4 


0.38+0.29 


-51+41 


- 


- 




RT.73 


-0.9 


-1.6 


6.94 


7.53 


12.0 


0.49+0.36 


-45+24 


- 


- 




RT.74 


57.2 


-64.0 


1.19 


1.72 


11.8 


- 


- 


- 


- 




RT.75 


31.0 


-9.1 


39.52 


48.67 


11.7 


0.63+0.12 


+64+7 


6.05+0.19 


-188+6 




RT.76 


28.3 


-6.4 


4.08 


5.32 


10.7 


- 


- 


- 


- 




RT.77 


25.2 


-61.4 


0.31 


0.69 


10.7 


- 


- 


- 


- 




RT.78 


24.0 


-60.4 


0.32 


0.67 


10.7 


- 


- 


- 


- 




RT.79 


28.0 


-27.4 


0.29 


0.41 


10.6 


- 


- 


- 


- 




RT.80 


-34.5 


8.4 


0.21 


0.64 


10.6 


- 


- 


- 


- 




RT.81 


26.8 


-41.0 


0.06 


0.08 


10.2 


- 


- 


- 


- 




RT.82 


29.3 


-1.6 


0.11 


0.13 


9.9 


- 


- 


- 


- 




RT.83 


29.9 


4.0 


0.45 


0.49 


9.7 


- 


- 


- 


- 




RT.84 


28.9 


19.1 


0.14 


0.14 


9.7 


- 


- 


- 


- 




RT.85 


21.7 


7.7 


2.72 


3.09 


9.6 


- 


- 


- 


- 




RT.86 


23.8 


2.1 


0.53 


0.88 


9.4 


- 


- 


- 


- 




RT.87 


26.2 


0.4 


0.42 


0.51 


9.4 


- 


- 


- 


- 




RT.88 


21.3 


3.1 


3.02 


3.57 


9.3 


- 


- 


- 


- 




RT.89 


21.8 


5.2 


0.64 


0.80 


9.2 


- 


- 


- 


- 




RT.90 


24.7 


3.8 


23.27 


28.03 


8.9 


0.11+0.01 


+38+3 


1.80+0.48 


-84 & +63* 




RT.91 


50.8 


-51.2 


0.52 


0.61 


8.8 


- 


- 


- 


- 


IRC+60370 


IRC.01 


-6.5 


14.5 


0.12 


0.13 


-39.7 


- 


- 


- 


- 




IRC.02 


-5.7 


16.1 


0.11 


0.12 


-40.0 


- 


- 


- 


- 




IRC.03 


33.2 


-35.7 


0.89 


1.10 


-40.2 


- 


- 


- 


- 




IRC.04 


-8.4 


13.7 


0.18 


0.29 


-40.5 


- 


- 


- 


- 




IRC.05 


4.7 


-10.0 


0.63 


0.75 


-40.6 


- 


- 


- 


- 




IRC.06 


-4.7 


15.8 


0.20 


0.24 


-40.8 


- 


- 


- 


- 




IRC.07 


14.3 


-10.3 


0.17 


0.22 


-40.9 


- 


- 


- 


- 




IRC.08 


12.4 


-9.9 


0.10 


0.16 


-41.3 


- 


- 


- 


- 




IRC.09 


-3.8 


16.3 


0.05 


0.07 


-41.3 


- 


- 


- 


- 




IRC. 10 


5.0 


-9.4 


0.16 


0.22 


-41.4 


- 


- 


- 


- 




IRC. 11 


17.2 


-8.8 


0.05 


0.06 


-43.6 


- 


- 


- 


- 




IRC. 12 


5.6 


13.2 


0.77 


1.00 


-43.7 


- 


- 


- 


- 




IRC. 13 


34.6 


-24.3 


0.08 


0.08 


-43.7 


- 


- 


- 


- 




IRC. 14 


28.1 


13.0 


0.08 


0.12 


-43.8 


- 


- 


- 


- 




IRC. 15 


-17.2 


-6.0 


0.11 


0.14 


-44.2 


- 


- 


- 


- 




IRC. 16 


18.6 


-10.5 


0.05 


0.05 


-44.2 


- 


- 


- 


- 




IRC. 17 


27.9 


12.3 


0.09 


0.16 


-44.4 


- 


- 


- 


- 




IRC. 18 


35.4 


-23.4 


0.11 


0.11 


-44.4 


- 


- 


- 


- 
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Table 7. continued. 



Source Feature 


"off 


Soff 


Peak Int 


Int Flux 


* peak 


Pl 


EVPA 


Pv 


B||[Gan.ij] 




(mas) 


(mas) 


(Jy/Beam) 


(Jy) 


(km/s) 


(%) 


O 


(xl0~ 3 ) 


(mG) 


IRC. 19 


6.5 


13.0 


1.03 


1.24 


-44.5 


- 


- 


- 


- 


IRC+60370 IRC.20 


3.9 


13.7 


0.10 


0.14 


-44.8 


- 


- 


- 


- 


IRC.21 


15.6 


-7.2 


0.17 


0.21 


-45.2 


- 


- 


- 


- 


IRC.22 


33.2 


-43.1 


0.06 


0.06 


-45.2 


- 


- 


- 


- 


IRC.23 


29.6 


9.9 


0.13 


0.21 


-45.3 


- 


- 


- 


- 


IRC.24 


2.8 


13.1 


0.08 


0.11 


-45.3 


- 


- 


- 


- 


IRC.25 


4.8 


11.0 


0.10 


0.14 


-45.6 


- 


- 


- 


- 


IRC.26 


4.0 


9.6 


0.11 


0.31 


-45.8 


- 


- 


- 


- 


IRC.27 


3.7 


8.4 


2.04 


3.29 


-47.0 


- 


- 


- 


- 


IRC.28 


4.5 


10.0 


0.05 


0.08 


-47.4 


- 


- 


- 


- 


IRC.29 


-7.6 


4.5 


0.15 


0.19 


-47.9 


- 


- 


- 


- 


IRC.30 


-8.9 


4.4 


0.20 


0.35 


-48.0 


- 


- 


- 


- 


IRC.31 


-0.8 


5.9 


0.10 


0.10 


-48.0 


- 


- 


- 


- 


IRC.32 


-5.6 


4.4 


0.08 


0.15 


-48.0 


- 


- 


- 


- 


IRC.33 


2.8 


8.0 


0.05 


0.06 


-48.0 


- 


- 


- 


- 


IRC.34 


27.8 


7.2 


1.14 


3.31 


-48.9 


1.51+0.10 


-132±2 


- 


- 


IRC.35 


28.3 


6.8 


6.29 


11.16 


-49.1 


0.61+0.07 


-67±7 


- 


- 


IRC.36 


28.3 


-5.7 


1.29 


1.47 


-49.3 


- 


- 


- 


- 


IRC.37 


3.6 


2.8 


1.19 


1.52 


-49.7 


- 


- 


- 


- 


IRC.38 


27.2 


6.5 


3.03 


4.35 


-49.8 


0.58±0.03 


-74±2 


- 


- 


IRC.39 


5.6 


5.8 


0.17 


0.24 


-49.9 


- 


- 


- 


- 


IRC.40 


30.6 


5.8 


0.22 


0.34 


-50.1 


- 


- 


- 


- 


IRC.41 


-1.1 


-0.4 


1.96 


2.54 


-50.3 


1.58±0.30 


-39±7 


- 


- 


IRC.42 


-0.9 


0.2 


11.91 


20.83 


-51.4 


0.65±0.30 


-58±16 


- 


- 


IRC.43 


4.4 


-8.6 


1.65 


2.05 


-51.8 


- 


- 


- 


- 


IRC.44 


0.0 


0.0 


51.23 


67.03 


-52.0 


0.57+0.02 


-77±2 


2.10+0.13 


+47±3 e 


IRC.45 


-0.6 


2.0 


5.21 


8.88 


-52.3 


0.45±0.07 


—93+5 


10.71 + 1.21 


+266±30 f 


IRC.46 


-1.3 


2.0 


1.68 


3.06 


-52.7 


- 


- 


- 


- 


IRC.47 


-2.0 


0.9 


1.65 


2.24 


-53.0 


- 


- 


15.52±4.05 


+331+82" 


IRC.48 


-0.8 


0.6 


11.26 


15.17 


-53.3 


0.46±0.18 


-75+14 


10.38+0.61 


+273±18 e 


IRC.49 


23.2 


12.6 


0.34 


0.48 


-53.8 


- 


- 


- 


- 


IRC.50 


23.3 


12.8 


0.33 


0.51 


-54.0 


- 


- 


- 


- 


IRC.51 


-1.6 


2.8 


9.08 


9.43 


-54.2 


0.19±0.02 


-97±6 


- 


- 


IRC.52 


-20.0 


13.0 


0.14 


0.26 


-54.6 


- 


- 


- 


- 


IRC.53 


16.7 


20.4 


0.06 


0.07 


-55.8 


- 


- 


- 


- 


IRC.54 


29.6 


3.9 


0.40 


0.49 


-55.9 


- 


- 


- 


- 


IRC.55 


29.2 


4.0 


0.09 


0.14 


-56.8 


- 


- 


- 


- 


IRC.56 


26.4 


21.2 


0.03 


0.04 


-57.4 


- 


- 


- 


- 


IRC.57 


28.8 


4.0 


0.23 


0.53 


-57.6 


- 


- 


- 


- 


IRC.58 


28.5 


4.6 


4.89 


5.48 


-58.3 


- 


- 


8.34±1.40 


-130±22 


IRC.59 


28.7 


3.7 


0.20 


0.29 


-59.0 


- 


- 


- 


- 


IRC.60 


14.6 


-4.8 


0.53 


0.57 


-59.5 


- 


- 


- 


- 


IRC.61 


14.4 


-4.5 


0.76 


0.81 


-60.6 


- 


- 


- 


- 


IRC.62 


7.0 


14.4 


0.07 


0.09 


-63.2 


- 


- 


- 


- 


e Edge/higher noise effects 


















* Blended feature 





















15 



